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. . . Abstract 
l> 

^>>^ ' In this work we consider superintegrable systems in the classical r-matrix method. 

0^ I By using other authomorphisms of the loop algebras we construct new superintegrable 

systems with rational potentials from geodesic motion on M^". 

>>. 

1 Introduction 

\o 

We shall consider classical integrable hamiltonian systems on the coadjoint orbits of finite- 
J> I dimensional Lie algebras according to [||. The dual space Q* to the Lie algebra Q is equipped 

' with the natural Lie-Poisson brackets specified by the condition that the Poisson bracket of 

two linear functions on Q* coincides with their Lie bracket in Q. Let -fT be a function on 
lO ■ 0*' ^ gradient of H. In the space C°°{Q*) of smooth function H determines the 

' evolution with the associated hamiltonian equation 

X = -{ad*^H) ■ X , xeQ*. (1.1) 

If Q is self-dual, i.e. has a nondegenerate inner product which allows to identify Q* with Q 
and ad* with ad, then (|1.1| ) takes on the usual form 

X = {H,x} = -{&d\7H) ■ X , xeQ. (1.2) 

Henceforth we shall always to identify Q* with Q and ad* with ad. 

Function I £ C°°(0) is called an integral of evolution with a hamiltonian H if 

{if,/} = 0. 

The evolution on a 2n-dimensional symplectic manifold M with the hamiltonian H is called 
completely integrable if there exists n functions Ji, . . . , which are independent integrals in 
the involution for the hamiltonian H 

{Ii,Ij} = 0, i,j<n. (1.3) 

The functions Ii, . . . , /„ are independent, if forms dli, . . . , din are linearly independent on the 
common level surfaces of these functions. 

The evolution on a manifold M, dimAf = 2n with the hamiltonian H is called superin- 
tegrable or degenerate, if there exists more than n independent integrals of motion {Ij}j^i , 
k > n and n of which are in the involution (|1.3[) M, O, p5[. For the superintegrable systems 



o 

C/3 



X 



all the integrals {Ij}j^i , k > n are generators of the polynomial associative algebra, whose 
defining relations are polynomials of certain order in generators (see [Q, |5| for a collection of 
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original papers). The main example of the superintegrable systems is a free motion with the 
following hamiltonian and equations of motion 



H = Y.Ph "ij = Pi ' Pj = ^^ (1-4) 

here Pj,qj are canonical variables on M = M^". Integrals of motion in the involution and 
additional integrals of motion may be defined as 

Ik = fiPi,---,Pn), Ijk = PjQk - QjPk ■ (1-5) 
Other known classical superintegrable systems with an arbitrary numbers of degrees of freedom 



are the harmonic oscillator, the Kepler problem and the Calogero system |21|. Notice, that 
the Kepler problem and the Calogero model may be obtained from the geodesic motion ( |l.4D 
on spaces of constant curvature Another examples of superintegrable systems can be 

constructed by using either purely algebraic techniques ^, ^, |l8| or separation of variables 
method at n = 2, 3 |2^, ^, ^. Some individual examples of superintegrable systems are listed 
in [21 1 with the corresponding references. 



Our aim is to show how superintegrable systems fit into a general pattern based on the 
notion of the classical r-matrix jl^, |2^. The main advantage of such embedding is that impor- 
tant structure elements for superintegrable systems, such as a Lax representation, separation 
of variables ^ and spectrum-generating algebra (dynamical algebra) |l^ can be system- 
atically derived from the underlying standard r-matrix formalism, which is a prerequisite for 
the study of the quantum case. 

We propose a dressing procedure allowing to construct the new superintegrable systems 
starting from known ones. As a natural initial point we shall select a geodesic motion ( |1.4[ ) 
on the Riemannian spaces of constant curvature. The Lax representations for these superin- 
tegrable geodesic motion are known |jl^ . To construct the new Lax equations associated to a 
potential superintegrable motion we apply the outer automorphism of the corresponding loop 
algebras |2|] directly to the Lax equations associated to a geodesic motion. 

The paper is organized as follows. In Section 2 we briefly recall notion of classical r- 
matrix method. By use a triangular decomposition of semi-simple Lie algebras, in Section 3, 
the algebraic approach to superintegrable systems is proposed and this scheme is applied to 
several examples. In Section 4 superintegrable systems are constructed in r-matrix formalism, 
while Section 5 contains some examples. 



2 Method of the classical r-matrix 

A systematic way for realizing integrable hamiltonian system on coadjoint orbits of the Lie 
algebras is provided by the r-matrix method |11, 22 1. 

Recall that the classical r-matrix on a Lie algebra is a linear operator R £ End(0) such 
that the bracket on Q 



[X,Y]j, = -i[RX,Y] + [X,RY]) 



X,YeQ, 



(2.1) 



satisfies the Jacobi identity |22, In this case there are two structures of a Lie algebra 
on the linear space Q given by original Lie bracket and by the r-bracket (^^, respectively. 
The Casimir functions tj on Q* invariant with respect to the original Lie structure are in the 
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involution with respect to the r-bracket. If r is an invariant function on Q*, the associated 
hamiltonian equation ( |1.1| ) on Q* is equal to 

^ = -ad*^.L, A=^R{dT{L)), L€Q*. (2.2) 
If Q is self-dual, then ( |2.2| ) takes on the usual Lax form |22|. 



It is obvious, that for any r-matrix in ( |2.2D all the Casimir functions give rise to integrals 
of motion in the involution |11, 22 1. We have to find the origin of an appearance of the 
special sup erintegr able hamiltonians and their additional integrals of motion. Application of 



the ad-invariant functions tj is a basic tool in the r-matrix method |11, 22| and, therefore, 
we consider these functions in greater detail to assume the standard identification of the dual 
spaces. 

To begin with let us recall some necessary facts from the notion of a universal enveloping 
algebra 0. Let be a Lie algebra and T{Q) be the tensor algebra of the vector space Q 

T = T^ QT"^ . .. , r'' = 9(g)9<g)...<g)0 n times . (2.3) 

If J be the two-sided ideal of T generated by the tensors 

x(g)y-y(g)x-[x,y], x,yeQ, 

then the associative algebra T/J is termed the universal enveloping algebra, which is usual 
denoted by U{Q). 

Let m > be an integer. The vector subspace of U{Q) generated by the products 
xiX2 • ■ ■ Xj, where xi,X2, ■ ■ ■ ,Xj G Q and j < m is denoted by Um{Q)- We have 

UoiQ) = C • 1 , C/i(0) = C • 1 © , U,{Q)Uj{Q) C U,+j{Q) . 

This sequence is termed the canonical filtration of U{Q). 

According by the Birkhoff-Witt theorem T{Q) = J(BS{Q) and algebra U{Q) is isomorphic 
to the symmetric algebra S{Q) as a vector space. If xi, 3:2, . . . , Xm G 0, then 

1 1 

w{xiX2 ■ ■ ■ Xm) = — X! P^^i^'i ' " ^-rn = — ^7r{l)a^7r(2) " " " X^r^rn) , (2-4) 
TTl. ^ TTl. ^ 

here Pt^ means the permutation operator corresponding to a certain Young diagram vr 
The map w ( ^.4| ) is a bijection of S{Q) onto U{Q), which is called the symmetrization. 

The Casimir functions or the ad-invariant functions on Q form a center 1(0) of U{Q). 
The symmetrization mapping ( ^.4| ) allows us to construct 1(0) by using the ad-invariants of 
commutative algebra S{Q) ||^]. 

The most interesting class of examples in the classical r-matrix method is provided by loop 
algebra C{Q, A) Algebra C{Q, A) can be realized as an algebra of the Laurent polynomials 
with coefficients in Q 

A0,A)=0[A,A-1] = |x(A)=^xAS xGsj 

and with the commutator [xX^,yX^ = [x,y]X'^~^^. 
According to an algebra homomorphism [l^] 

U{CiQ,X))^C[X,X-']0UiQ) 
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the Casimir functions on the loop algebras can be recovered by the ad-invariants Tj{x) in 1(0) 
r,-^ = ResU^o'A(A)-r,(x(A)), r,-(x)Gl(0), </)(A) G C[A, A-^] , (2.5) 



where (p{X) is some rational function on spectral parameter A with numerical values [ p2[ . 
Studying the superintegrable systems we have to introduce the usual tensor algebra 



r-(0, A, i/) 



£(0,A)®£(0,/i) 



m — times , (2.6) 



and canonical filtration of the corresponding enveloping algebra U{Q, X, fi, . . .) generated by 
subspaces Um{Q, A, /i, . . . , i/). These vector subspaces are produced by subspaces UmiQ) 

•^i\i2---ik • • • ! '^) — ^ ^ Xi^i2---if,X'^ fl'^ ■ ■ ■ U'^'' , k m , Xjj^jj - ife ^ ^m(0) • 

jl,j2,---,jk 

In just the same way as for a Lie algebra 0, we can define the canonical mapping of the 
loop algebra C{Q, A) into U{Q, X, fi, . . .). Any element L{X) of C{Q, A) can be embedded into 
Um{Q,X,l^,---,'^) 



Lj{Xj) = icii © • • • © © L{Xj) © idj+i © • • • © idm G Um{Q, X, fi, . . . ,iy) . 



and 



^2, . . . , A,) = n ^.n(AjJ , 1 < A; < 



n=l 



L^"^) (A, ^, . . . , I/) = L(A) © L(/i) © • • • © L(i/) 



m — times . 



(2.7) 

(2.8) 
(2.9) 



here Ai = A , A2 = ^, ... Xm = 

If Q is identified with its dual, then r-bracket (2J.) can be rewritten in the tensor form 
|ri| , 22, Let L{X) G X) be a generic point in the loop algebra, which is regarded as a 
Lax matrix. The corresponding r-bracket is given by 



{Li(A),L2(^)} = [ri2(A,M),^i(A)]-[r2i(A,/.),L2(^)], 



7^21 (A, ^) = A2n2(A,/U)Pi2 



(2.10) 



where P12 is a permutation operator in £(0,A) © £(0,/i) and rij{X,fi) are kernels of the 
corresponding operators R and R* in ( |2.1D [p^]. Notice, that the r-matrix scheme is extended 
easily to the twisted subalgebras of loop algebra C{Q, X) and the corresponding matrices ri2 
have rational, trigonometric and elliptic dependence on spectral parameter. 

The Poisson brackets between the elements ^jlj2 - jk^^^ /^j • • • > i^) ( |2.8| ) can be written in 
the "generalized" r-matrix form. For instance 

{Li2(A,//),L3(i/)} = [ri3(A,i/) + r23(^,z^),-Li2(A,/u)] 



(2.111 



[r3i(A,z^),L23(/i,i^)] - [r23(/U,i/),Li3(A,z^)] 



where rjj(Aj,Aj) are r-matrices, which act nontrivially in the corresponding subspaces of 
Tm{Q,X,iJ,,...,u) (l^l) and 



{Li2(A,/i),L34(i/, ry)} 



rW(A,//,i/,r/),Li(A)] + [r(2)(A, ^, z., ry), L2(/i) 
r(3)(A,^,i/,r?),L3(zy)l - [r(^)(A, ^, i/, r/), ^4(7?) 



(2.12) 
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where 



(Dm 


rj) 




-P34 , 




ri) 


= Pi2r«(^,A,i.,?7)Pi2, 




r(3) (A, 


rj) 


= Pi3 IH^, i^, v) + Pi2r{^i, A, 


??)Pl2] -Pis 




rj) 


= P34r-(3)(A,/i,r?,i.)P34, 






rj) 


= n3(A,J^)-L24(/i,??) • 





Here Pjj are operators of pairwise permutations in the tensor algebra Tm{Q, A, /x, . . . , z^) 

Integrals of motion in the involution are completely defined by the ad-invariants ( |2.5| ) 
of £(0,A) and by the linear r-bracket ( |2.1[l| ) [ p] , p^ . For description of the superintegrable 
systems we have to consider ad- invariants of f7(0, A, /i, . . .) and more complicated embedding 



"^7iL- -7 (-^'/^'••- (2.8-2!9|). In the next section we consider several examples, that allows 



us to understand the origin of the appearance of additional integrals of motion in the classical 
r-matrix method. 

3 One class of superintegrable systems 

Let be a simple Lie algebra with a Cartan subalgebra b and a system of simple roots 
n = {ai, 02, • • • , On}- Let 

— S-oti ) — 1 fi — ^—cii ) ^ — 1 . . . Tl , 

be the Chevalley generators and {hi, Cj, fi\ be a Cartan- Weil basis @, normalized by (e^, /j) = 
1: 

[hi, hj] = , [ei, fj] = 6ijhi , 

[hi,ej] = CijBj , i^iT fj] ~ ^ijfj ; (^'1) 

(ade,)-^'^+i -6^ = 0, (ad/i)-^'^+i ■ fj = . 

Here Cjj are entries of the Cartan matrix, which are integers in this normalization. 
Let n+ (tT_) be the linear span of root vectors ej {fj) such that 

g = n_eben+. (3.2) 

This decomposition is termed the triangular decomposition of Q. 

Every basis in the commutative subalgebra Um(b) C Um{Q) could be associated to the 
family of functionally independent integrals in the involution. Let us define a hamiltonian H 
as a function on alone generator hi of the Cartan subalgebra b 

H = H{h,) , /i, G b , H£ U^{b) . (3.3) 

It is immediately seen that the n — 1 functionally independent integrals in the involution 
may be constructed for the hamiltonian H (|3.3D in Um{b) 

Ij=Ij{hi,...,hn), hi eh, Ij£Um{b), (3.4) 

here n = dimb = rankQ. 
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Moreover, we can introduce another n{n — 1) integrals of evolution in [7(n_ © Tl-i-) 



4^ = e^.-/rGC/p(n_©n+), ^ = ^, p = max{/,m}, 



(3.5) 

dH 



where I > 1 and m > 1 in monomials /j™ are positive integers in normalization (js] 

Let p > 1 be a largest power of monomials /j^ , j,k = 1, . . . , n. Taking ( |3.1| ) into account 
we observe that other monomials in Up{Q) 

T^. = h^^egfg, n+j2 + j3<P, (3.6) 
belong to the generalized dynamical algebra |5| defined as 

dH 



{H{h),Tllill] 



In this case the generalized dynamical algebra ^ of (i]l) generated by elements (| 
coincides with the polynomial subalgebra Up{Q). Operators 



\dhiJ 

are counterparts of the ladder operators of hamiltonian H ( |3.3|) in quantum mechanics. If 
the hamiltonian H is a linear function on generators, the spectrum-generating algebra 19] 
generated by elements (^) coincides with the polynomial subalgebra Up{Q). 

The realization of a proposed scheme for a given hamiltonian is a quite difficult task, 
which is not always dictated by an a priori obvious procedure of the searching of a necessary 
algebra Q. Moreover, the different algebras Q could be associated to the single superintegrable 
system. 

On the other hand, taking some fixed algebra Q, we have to introduce a suitable repre- 
sentation for a realization of superintegrable system in canonical variables on it's coadjoint 
orbits. 

Let us take the n three-dimensional Heisenberg algebras f)^- {Q = ©f)j) having basis 
{hj,ej,fj) such that [ej,fj] = hj and [hj,ej] = [hj,fj] = 0. The center of nilpotent algebra 
f)j is C ■ hj. The corresponding universal enveloping algebra C/(0 = ©f)j) is a Weil algebra 
Wn- Polynomial subalgebra U2{Q) of Wn is well-known in mathematics [||, 0] and physics 
IQ, ^ |T^. It consists of quadratic generators Tij, Pij, Sij and linear generators {hj,ej,fj) 

Tij — Cj/j , Sij — 6j6j , Pij — fifji (^-9) 

where we do not consider polynomials in elements of the center of f)^. 
The non- vanishing polynomial Poisson brackets in U2{Q) are 

{/i; — hiSijfm i {fit Sjm\ — hi {^^ijdi-n ~\~ ^im^^j) ? 

{Tij , Tim} = i^imhiTij — SjlhjTim) , (3.10) 

{Tij, Sim} — ~hj {PjmSil + SjlSim) i {Tij, Pirn} — hi {dimPjl + ^ilPjm) , 

{Sij,Plm} = hj {SjmTil + SjlTim) + hi {6imTjl + SuTjm) , 
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The coadjoint orbits of the Heisenberg algebras are planes hj = aj £ C. On the fixed orbit 
of the algebra Q = ©f)^ we can reinterpret the brackets ( p.lOl ) not as polynomial brackets 
in generators Q, but rather as linear brackets between generators T^j, Pij, Sij and ej, fj 
with the structure constants hj = aj. It allows us to identify representations of U2{Q) with 
representation of the known Lie algebras. If we consider n equivalent orbits of nilpotent 
algebra f) (oj = aj), the complete algebra U2{Q) is isomorphic to the semi-direct product 
Wn sp{2n) 1^, ^. Subalgebra of generators Tij is isomorphic to the algebra u{n) or, after 
separation of Casimir operator J2'^kk, to the algebra su[n). 

The n-dimensional commutative ideal in U2{Q) (the Cartan subalgebra) is generated by 
elements hj = ejfj = Tjj. If we introduce the hamiltonian H by the rule ( |3.3| ) on this 
n-dimensional ideal, then another generators of U2{Q) have the following brackets with it 

H = J2i=i hj = YJj=i Tjj , {H, Tij} = {hi - hj)Tij , 

{H, ej) = -hjCj , {H, fj} = hjfj , (3.11) 

{H, S,,} = -{hi + hj)Sij , {H, Pi,} = {h, + hj)P,j . 

The hamiltonian H is a sup erintegr able hamiltonian in f/2(0) on the following coadjoint orbit 
of = ©()j- 

hi = hj, i,j = l,...,n. (3.12) 
Note, that we can select different n-dimensional commutative ideals in the various subalgebras 



Um{Q) to construct superintegrable systems (3^3^) with some fixed hamiltonian (3.3). 



The n-dimensional harmonic oscillator in the euclidean coordinates is described by the 
hamiltonian 



The phase space of oscillator ( p. 13 ) can be identified with the coadjoint orbits of Heisenberg 
algebras according to 

Bj = aj = ^{pj - iajQj) , 

fj = a+ = -^{pj + iajQj) , (3.14) 
hj = iaj , 

where aj, are the standard creation and annihilation operators. In this realization, on 
the orbit ( |3.12D hamiltonian H (3.13) of isotropic oscillator is superintegrable in U2{Q) 0]. 



However, it will be useful to identify the phase space of oscillator with the coadjoint 
orbits of another algebras. For instance, if the phase space of oscillator be identified with the 
coadjoint orbits in sl{2)* , we can introduce integrable systems isomorphic to oscillator 

H' = \T.{p] + -hf)+Y.'f2^ (3.15) 

^ 3=1 i=l 



which is known as the Smorodinsky-Winternitz system at n = 2, 3 |2c 

Let us prove it. Consider an infinite-dimensional representation W of the Lie algebra 
sl{2) in linear space V defined in the Cartan- Weil basis {/i, e, /} in End(V) equipped with the 
natural bracket 

[h,e\=e, [h,f] = -f, [e,f]=2h, A = h' + ^{ef + fe) . (3.16) 
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If operator e is invertible in End(V), then the mapping 



h' = h, 



(3.17) 



is an outer automorphism of the space of infinite-dimensional representations of sl{2) in V 
|26|. The mapping (3.17) shifts spectrum of A on the parameter /3 

A ^ A' = A + /?. 

We can suppose that automorphisms ( 3.17| ) defines one-parametric reahzation W{P) of s/(2). 
For instance, reahzation of s/(2) with one free parameter /3 in the classical mechanics is given 
by 



2 ' 



(3.18) 



where {q,p) is a pair of canonical coordinate and momenta with the classical Poisson bracket 
{PtQ} = 1- Motivated by realization (3.18) we present one application of the automorphism 
(|3.17] ) in the theory of integrable systems. Consider a classical hamiltonian system completely 



integrable on the M " with the natural hamiltonian 

n 



,qn) 



Let the phase space be identified completely or partially with the m coadjoint orbits in s/(2)* 
as ( P^ . Then the mappmg 

A- 



H ^ H' = H + }_^^ , Pj eR (3.19) 

j=i 

preserves the properties of integrability and separability. The list of such systems can be 
found in |21]. 

Let us take the n algebras sl{2) {Q = ©s/(2)) having basis {hj,ej,fj) with the brackets 



[hj,ej] 



[hj, fj] = -ajfj , [ej, fj] = 2ajhj , 



(3.20) 



where {aj}'j^i G C is a set of arbitrary constants. 

As usual, we can define the polynomial subalgebra U2{Q) of the corresponding universal 
enveloping algebra. Here we consider a special subalgebra of U2{Q) consisting of quadratic 
operators Tij G U2{Xl- © n+) , i ^ j and linear operators Tu G Ui{b) 



hi , 



(3.21) 



These operators are closed under the Poisson bracket 



or, in the unit polylinear form. 



2 iyOijbjjYl 

— OLkSik) TjmTik , j ^ k, I ^ m , 
Oj {6jl - 6jm) Tim , I ^m, 
0, 



(3.22) 



{Tjk,Tim} — [oijSjk {Sjl — Sjm) Tim + OllSml {^kl — ^jl) Tjk + 

+ 2 (1 - 6jk){l - Sim) {aj6jm - akSik) Tj^Tik ,] 
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As before, we introduce linear hamiltonian H = J^^j — J2 '^jj Cartan subalgebra of 

Q = (Bsl{2). Elements Tij of U2{Xl^ ©^1+) have the following brackets with the hamiltonian 

{H,Tij} = (oj - ai)Tij . 



The coadjoint orbits in sl{2) algebras ( 3.20| ) are fixed by the value of the Casimir operator 



A ( 3.16| ): Aj = a| G C. On the special coadjoint orbit of = (Bsl{2) determined by 



Oi = aj , i, j = 1, . . . ,n the linear hamiltonian H (|3.3|) has additional integrals of evolution 



The phase space of oscillator can be identified with the coadjoint orbit of Q by using the 



representation ( 3.18 ) and the similar transformation 



Matrix elements {hj,e'j,fj} in ( ^.23 ) are generators of sl{2) in representation ( 3.1S| ). More 
explicitly 



1-2 „2„2 , r,„-„ „„\ ^^„+\2 



ej = -{pj - ajQj + 2iajPjqj) = ^{ap , 
fj = \ip] - - '^ioLjPjqj) = ^ «j , 

hj = :^{p] + a]q]) = ^N,, (3.24) 



where {aj, a^' , Nj}'^^^ are the standard operators ( p. 141). H amiltonian H = 2iJ2 Tjj of oscilla- 
tor ( |3.13| ) or of the Smorodinsky-Winternitz system ( |3.15| ) and integrals of motion Tij ( ^.21 ) 
are at most second order polynomials in generators of sl{2)* . 

If the ratio of frequencies aj/a^ = l/ra is rational, then the hamiltonian of oscillator 
( 3.13| ) has additional integrals of motion in Up{Q) (|3.5D |||]. The corresponding dynamical 



algebra coincides with Up{Q). For instance, the anisotropic oscillator in two dimensions with 
a 2 : 1 ratio of the frequencies and with a third order polynomial algebra U^{Q) has been 
considered in p. 

So, we can see that for the superintegrable hamiltonian H ( |3.3D in Um(^) the additional 
integrals of motion belong to Um{^- © TI+)- In the next section, motivated by the presented 
examples, we consider superintegrable systems in the classical r-matrix method. 

4 Superintegrable systems on sl{n) 

In application to integrable systems every element L{\) of the loop algebra £(0, A) be a matrix 
in some fixed matrix representation of the algebra Q and representation space be an auxiliary 



space [11|. Let us consider an algebra Q = sl{n,C) in the fundamental representation and 
begin with the standard rational r-matrix ||ll|, |2^. In this case element L{X) of C{Q,X) be 
a n X n matrix in the auxiliary space C" and the matrix elements of L(X) are some rational 
functions of spectral parameter A. 

Basis of invariant functions in I{sl{n, C)) could be selected as 

Tk{L{X)) = hvL^X) , k<n. (4.1) 

The family of the integrals of motion in the involution is generated by Tk{L{X)) ||Tl|, ^ 

/,,fc(L) = cD«(rfc(A)), (4.2) 
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where ^>1^*^ are various linear functionals defining a set of the functionahy independent integrals 
of evolution. For instance, 

= Res|^^o(</.,(A)-z), 0(A) eC[A,A-1], (4.3) 

here 0i(A) are some functions of spectral parameter with numerical values. In this case 
integrals are at most /c-order polynomials in generators Q, li^k £ Uk{Q)- 

The Lax equation (p.2|) associated with the hamiltonian Jj^^ ([4.2|) is equal to 



dt 

where the second Lax matrix Aij. has the form 



"^^^^^ -- {hk,m} = [L(/.),^,„,(/.)] , (4.4) 



Ai^kif^) = cD^tn (r2i(A, /.)L^i(A)) , (4.5) 



here trace tr^ is taken over the first auxiliary space. 

The representation space of the subalgebra Um{Q, A, /i, . . . , i/) be an extended auxiliary 
space 

y{m) ^ ^nr^^y, = Vi ^ V2 ^ ■ ■ ■ , ~ C" . (4.6) 

The elements Lj (Xj) ( |2.7| ) and L^'lj^-'-jk ■ ■ ■ ^ ^) ( 1^-^ ) belonging to the algebra C/m(0, • • • ) ' 
are the x matrices in y^™) 

Lj{\j) = /i • • • L{\j) (g) Ij+i (g) ■ ■ ■ (g) Ijn , (4.7) 

m 

L^™) (A, ,1/) = J]^ Lj(Aj) , Ai = A, A2 = A*, . . . , Am = z/, 

i=i 

here / means a nxn unit matrix and subscript j shows in which of the spaces Vj in the whole 
space y^*") the matrix L{X) acts nontrivially. 

The equation of evolution for the matrix L^"^) has a commutator Lax form 

(A, /...., I.) = [l(-) ( a, ^ . . . , I.) , A(-) (A, M . . . , ^^)] , (4.8) 
with the following second matrix 

m 

A^'^\X,H ...,i,) = Y,Aj{X,), Ai = A, X2 = ^l,...,Xm = l^ (4.9) 

which is a sum of the matrices Aj{Xj) of the type ( |4.7| ) acting in the whole spaces 

The spectral invariants of the matrices L^™) give rise to an involute family of the integrals of 

motion as before. 

Taking into account the symmetrization mapping w (^^) let us introduce matrices 

l("'^)(A, fi,...,i^) = P.l(™)(A, fi,...,i^) = P.Li(A)L2(//) • • • LM . (4.10) 
Here permutation matrix P^^ in V^"^^ is determined by 

(xi X2 • • • Xm) = g> x^^2) ® • • • ® > (4-11) 

for any set of vectors xj in or 

P,-AiP2---^™ = ^.(i)i?.(2)---^.M-P., P' = ^, (4.12) 
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for any n x n matrices A, B, . . . , D embedding in according to (4.7). The permutation 

of subscripts in ( [4.11 ) and ( 4.12| ) is defined by a certain Young diagram vr |15]. 
The equations of evolution for the matrices L^™''^) are equal to 



dt 



L^'™'") (A, . . . , z.) = (A, ^ . . . , I/) - (A, 



,iy)L 



(m,7r) 



(4.13) 



where matrix A^^^ is given by (|4.9| ) and the second matrix A^"^''^^ differs from it by the 
permutation of spectral parameters in accordance with the Young diagram vr 

(A, ^ ...,!/) = P^^('"'-)(A,^ ...,z.)p-\ 

m 

= X! (^tO) ) ' ^1 = ; A2 = /i, . . . , Am = 



The right-hand side of ( 4.13 ) is a matrix commutator if and only if 

A(-)(A,/....,z.) = A(-'-)(A,M...,i^), ^ ^(A,) = ^(A^o-)). 



(4.14) 



Assuming ( [4. 14]) holds we can define new multivariable generating functions of the integrals 
of motion 



m ^ ' 



(4.15) 



^s^(A,/i ...,i/) =0, 



The trace tr^^) in ( 4.15 ) is taken over the whole space 

Vim), 

For any hamiltonians li^k ( [4. 2D condition (4.14) is always fulfilled for the equivalent spec- 
tral parameters 

A = Aj = A^(j) , j = l,...,n, Vvr (4.16) 

that corresponds to a choice of another basis of ad- invariant functions in the center l{g) ||l^ 
by using the outer powers of matrix L{X) 



(A) = ltr(„)L('"'-)(A,A,...,A). 



(4.17) 



Symmetric functions Sm(A) can be expressed in the symmetric functions Tm{X) (|]^) according 
to Newton's formulas |]15[ and functions Sm{X) give rise to integrals in the involution as before. 

In additional, condition ( [4.14| ) is fulfilled for the independent on spectral parameter matrix 
A in 

AiM = ^?tri (r2i(A,/u)L™-^(A)) = const / 0. 

We observe that matrix Ai^mifJ') depends on spectral parameter // via r-matrix only. Hence, 
the constant in spectral sense matrix A(^) and the corresponding special hamiltonian are 
related to the singular points of the r-matrix. 

The standard rational r-matrix for algebra sl{n, C) is equal to 



'^12(A,^) 



Pi 



12 



A — 

We can choose the special linear functional in ( [4.5| ) as residue at infinity 

^x{z) = - Resl;,^^ (</.(A) • z) , 



(4.18) 



(4.19) 
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such that the second matrix 



A = tri 



P- hm (/)(A)L™-XA) 

A^oo 



(4.20) 



is independent on spectral parameter. Moreover, if some invariant polynomial rfc(A) ( f4.1[) has 
a nontrivial residue at A = oo on the phase space 



H 



Res I 



A=oo 



(4.21) 



which is chosen as a hamiltonian, the corresponding Lax matrix A ( ^.201 ) do not equal to zero. 
In this case the singular point A = oo of r-matrix ( 4.18| ) is associated with the superintegrable 
hamiltonian H ( 4.21| ) and with the new multivariable generating functions of the integrals of 
motion ( 4.15| ) 



s-(A,^,...,i/) = ltr(^)L(-'-)(A,/i 



(4.22) 



We may to lose the property of involution for these new multivariable functions and the 
corresponding integrals of motion 



{sJ^"(Ai,A2,...,Am),Sfc''(Aii,/U2,---,Atfc} + 0. 



(4.23) 



These brackets are completely recovered by the polynomial r-brackets for the matrices L^^''^^ , 
as an example see (|2.11[ - ^?T^ ). 

The coefficients of the spectral curve of L{X) determined by the characteristic equation 



C{z,X) = det {zI + L{X)) = 0. 



(4.24) 



give rise to integrals in the involution only. Complete set of integrals is determined by the 
generalized spectral surfaces 



C{z, ... u) = det (zl + (A, /i ... I/; 

As an example, let us consider the 2x2 Lax matrix 



0, 



m < n . 



L(A) 



h e 
f -h 



(A), 



(4.25) 



(4.26) 



in the two-dimensional auxiliary space with a suitable basis |11] 



P 



1 \ 

10 

10 

\ 1 / 



where Pauli matrices aj form an orthonormal basis of ■sl{2) in C^. 
Matrix L{X) has alone invariant polynomial ( 4.17| ) 

S2(A) = ^tr [PL{X) (g) L{X)] = detL(A) = -^tiL'^{X) 
and one second order complementary polynomial ( 4.22| ) 



:T2(A), 



.,(A, ,) = Itr [PLiX) ^ m] = HXM,) + <A)/(,) + /(A)e(^) 



(4.27) 



(4.28) 
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These polynomials may be the generating functions of the integrals of motion for some super- 
integrable hamiltonian (compare with ( |3.5D ). 

Next let us consider the Belavin r-matrices 0, which are meromorphic solutions of the 
classical Yang-Baxter equation such that 

n2(A, fi) = ri2(A -fi) = -r2i(/i - A) , (4.29) 

Pl2 

ri2(A-/i) = - hO(l), Pi2X®y = y®x. 

A — H 

Here r^j are rational, trigonometric or elliptic matrix- function on spectral parameters [Q, 11, 
22|. These matrices have a pole at A = ^ with a permutation operator P12 ( [1.29| ) as residue. 
In additional the rational function r(A — fx) of two variables A and /i has the special point 
at A = 00 in its domain of definition. In this point function r(A — /i) has a distinct from 
zero residue, which is independent from second spectral parameter /x. The generally accepted 
elliptic r-matrices |ll|, |2^] have not such special points in its domains. Nevertheless, the 
similar to the rational case construction can be proposed as well. 

The Lax equation (2^) and r-bracket ( |2.1[l| ) are covariant under the similar transforma- 
tion 

L U-^LU , A U-^AU , nj U^^U2^rijUiU2 , (4.30) 

where C/ is a constant matrix on a phase space. If the matrix U depends on spectral parameter, 
we can use such similar transformations to construct the additional poles of r-matrix. In this 
case the multivariable generating functions s'!^{X, fi, . . . ,1^) are changed 

•Sm(A, /i, . . . , I/) ^ sl,{X, fi,...,iy) = — tr(™) [Zt,L{X) ^ L{fi) ...(g) L{u)] , (4.31) 

here projector in ( [4.22| ) is substituted by matrix 

= [U{\) ® U{fi) ... U{u)r^ ■ ■ U{\) ® U{fi) ... U{u) . (4.32) 

depending on spectral parameters. For the equivalent spectral parameters functions ^^(A) 
are covariant under the similar transformations. 

As an example, here we study an elliptic r-matrix on the twisted loop algebra £(s/(2), a). 
Let us consider the period lattice T = 2K7L + 2iK'1j, where K and K' are the standard 
elliptic integrals of the module k £ [0,1], and introduce the corresponding elliptic theta 
function &ij{X,k) as in [Q. In these notations the standard elliptic r-matrix is 

3 

r{X- ij.) = ^Wk{X- ■ aki^ (Jk, (4.33) 

k=l 

where are the Pauli matrices and 

^ 9;i(0,A:)9io(A,fc) ^ e[,{0,k)eoo{X,k) ^ e[,{0,k)Qoi{X,k) 

"^'^ > eio(0,A;)en(A,A:) ' ""'^^ 600(0, A;)eii(A, A:) ' '^^^> 601(0, A;)en(A, A;) ' 

Function r(A — fi) ( |4.33D is meromorphic in C and has simple poles at A = ;U modF. 
According to we introduce the similar r-matrix 

p(A, fi) = C/i2^(A, fi, Aoo) r{X - fi) Ui2{X, fi, X^o) , 
Ui2{X, fi, Aoo) = f/i(A - Aoo + K)U2{fi -Xoo + K). 
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with the following matrix U{X) 
UiX) = 



Goi e 

-9oo —©01 



(4.34) 



here Aoo is an arbitrary point. This r-matrix p{X, jj) has been introduced for the purpose of 



separation variables in [24|. More explicitly 



i=i 

+ [wsifJ, - Aoo, A;) - ws{X - Aoo, A;) +wo]-(7s0as 



(4.35) 



+ iw{lJ. - Aoo, fc) • (T3 (g) (72 - iw{X - Aoo, ^) • 172 ® 0-3 ; 



w(A, k) = wi{X, k) = W2{X, k) 



©;i ©io(A,fc) 

910 Gil (A, A:) ' 



u;3(A) 



e;i(A,fc) 

eii(A,A:) ' 



here wq is a some constant |24]. The poles of p{X,p) by the first spectral parameter A be at 
the points 



A = /i mod r , 
A = Aoo niod r , 



Res|;,^^/)(A,^) = P, 

Resl;^^;^^ p{X, ^) = Z = (0-3 + cj_ + cr+) (g) 0-3 . 



So, if L{X) is an orbit of the new r-matrix p(A,/x), such that the second matrix 



A = tri 



Z- lim (/'(A)L(A) 

A *Aoo 



is nontrivial constant in spectral sense, the special hamiltonian 



H = Resl tr 



mL\x) 



is sup erintegr able. 

In the next section we present some nontrivial examples of r-matrix orbits associated with 
superintegrable systems for the rational r-matrix. 



5 Examples 

The above construction of superintegrable systems can be applied to the Gaudin magnet, 
which was introduced in the quantum mechanics The classical version turned out to be 



useful example for developing a general group-theoretic approach to integrable system |11, 22]. 

We shall consider the rational Gaudin magnet related to sl{N) algebra. The model in 
question is defined on the M coadjoint orbits of sl{N)* in variables X^™'\ (m = 1, 
i,j = 1, . . . , A^). The corresponding Lie-Poisson brackets are 



The coadjoint orbits are fixed by values t^™^ of the ad-invariant functions on sl{N) 



(5.1) 



(5.2) 
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The Poisson bracket (^]^) is nondegenerate on the manifold ( |5.2D having dimension n = 
MN{N — l)/2 for the case of generic orbit (all t'f^^ are distinct). In what follows we assume 
that the orbit is generic. 

Fixing some element Z G sl{N) as a residue at infinity we consider the special Lax matrix 
L(A) E C{esl{N)) 

L{X) = Z+Y.^—^, (5.3) 

m=l 

where {5m} is a set of M arbitrary constants. Matrix L{X) obeys the linear r-bracket ( p. 10 ) 
1 11, ^] with the rational r-matrix (|4.18| ). 

The basis elements t/. (4.1) are meromorphic functions of A 

M k jj 

rfc(A) =Ck+Y.Yl (X )j ' 

m=lj=l ~ 

here = k^^iiZ^ and I^k ~ ^'^^ fixed constants. Another residues form a family 
of n = MN{N — l)/2 independent integrals in the involution. It is immediately seen that the 
special hamiltonians in this family 

M 
m=l 

are nondegenerate functions on the generic coadjoint orbits of sl{N)* and they are corre- 
sponded to the constant in spectral sense second Lax matrix 



j^{k) _ ^k-1 ^ 

So, hamiltonians H^^'^ are superintegrable hamiltonians and complete set of the integrals of 
motion can be generated by ( 4.22 ). As an example, additional integrals of evolution may be 
constructed from the quantities 

Cfc = Res s;;(A,^,...,i/). (5.4) 

Here Res means residue of some fixed order k = (ki, k2, ■ ■ ■ , km), kj < Kj at the points 

X = 5j^, fi = , . . . , u = Sj^ , 

which belong to the divisor of the poles of multivariable function sJ^(A, fj., . . . ,1^) 

D = {{6j,K,),j = l,...,M,{oo,l)}. 

As a second example, let us consider superintegrable natural systems on M^"" with the 
following hamiltonian 

H = T +V = '^aijPiPj + V{qi,. . . ,qn) , UijGR, (5.5) 



where {pj,qj}'j^i are canonical variables. By V = in (|5.5| ) this hamiltonian of geodesic 
motion is superintegrable ( |1.4| ). In the classical mechanics the geodesic motion on the Rie- 
mannian spaces of constant curvature is isomorphic to the rational Gaudin magnet on the 
algebra Q = ®s/(2) [l^. Taking the known Lax matrix L{X) in C{sl{2), X) associated to a free 
geodesic motion we can get an infinite set of the Lax matrices associated to the completely 
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integrable potential systems |26]. Construction of the new Lax equations consists of the appli- 
cation of the outer automorphism of infinite-dimensional representations of sl{2) ( |3.17 ) to the 
loop algebras. As a second step, we have to apply the consequent projections of this mapping 
onto a suitable Poisson subspaces of r-bracket (ad|j-invariant subspaces). The superintegrable 
systems are related to the special projections, which lead to a constant in spectral sense Lax 
matrix ^(A) (U). 

Let us consider the Lax matrix L{X) in £(s/(2),A) associated to geodesic motion. Here- 
after we propose that this matrix L{X) in the form ( |4.26D obeys the r-matrix bracket ( p. 10 ) 
with the rational r-matrix ( 4.181) . Following to the general procedure js], [ij, 26] let us intro- 
duce new Lax pairs 



L'(A) = L(A) -(T_ • U(A)e-^(A) 



MN 



(5.6) 



^'(A) =A-a-- ^{X)e-\X) 



MN \ UmnW 



Here 0(A) is a function on spectral parameter and [2:] mat means restriction of z onto the ad|j- 
invariant Poisson subspace of the initial r-bracket |14, For the rational r-matrix ( [4.18 ) 
we can use the linear combinations of the following Laurent projections 



\z\mn 



+00 



N 

E 



ZfcA 



(5.7) 



MN 



'M 



or the Taylor projection by M = 0. 

The mappings (|5.6|) play the role of a dressing procedure allowing to construct the Lax 



matrices L^,^^(A) for an infinite set of new integrable systems starting from the single known 



Lax matrix L{X) associated to one integrable model. The Lax matrix L'{X) (5^) obeys the 



linear r-bracket ( p. 10 ), where constant 
dynamical variables 



-matrices substituted by r^ -matrices depending on 



{[4>{X)e-\X)\MN - [<i}{n)e~\lj)]MN 
(A-^) 



• (7_ 



(5.8) 



The Poisson map 



\>2n 



equations of motion |14]. Letting 



£(0) may be determined by using a Newton form of the geodesic 



A 




(5.9) 



which is independent on spectral parameter A, we observe that the Lax matrix L{X) ( [4.26| ) is 
recovered by the alone entry e(A) and by the hamiltonian H according to 



L(A) 



(A), 



(5.10) 



In this case the generating functions S2{X) and S2(A) obey the following equations of motion 



ds2{X) 
dt 

ds'^jX) 
dt 



0, 



0, 



dleiX) 



0, 



1 o 1 

-dx + umnWOx + ^umn,xW 



(5.11) 

e(A) = SiK/iv(A)]-e(A) = 0. 
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Here Bi[umnW] is the Hamiltonian pencil operator for the coupled KdV equation 



superintegrable systems correspond to the special solution of ( 5.11 ) by OumnW/QX 
integrate form of the KdV recursion relations ( ^.11 ) are 



|. The 
0. The 



4(A) 



e • e. 



e • e 



(5.12) 



A simple substitution for the entries of matrix L(A) 

e(A) = B\ h{X) = -ej2 = 
fix) = -e,,/2 = -Bl-BB,,, 



-BB, 



(5.13) 



turns determinants ( ^.12]) into the form 

S2{X)=B^B^,, s'^iX) = B^B^^ + B^ 



^4 



MN 



if we use an explicit formula for the potential UMNiX). These equations have the 
Newton's equations for the function B 



B.. 



B,. 



S2{X)B- 



4(A)S" 



(5.14) 
form of 

(5.15) 



B 



/(A) 



MN 



If we assume that B = J2jLo QN~j X^ is a polynomial, then its coefficients qj obey the Newton 
equation of motion ( 5.15| ) with S2(A) = J^^kX^, where Ik are integrals of motion. Here we 
reinterpret the coefficients of S2{X) and S2(A) in ( 5.15| ) not as functions on the phase space, 
but rather as integration constants. In variables qj mapping (^^) affects only on the potential 
(g-dependent) part of the integrals of motion J^. The kinetic (momentum dependent) part 



of Ik remains unchanged. So, the dressing mapping ( ^.61) allows us to get over from a free 
motion on M^" to a potential motion on M^". 

To construct the multipole Lax equation we introduce an appropriate completion Cd{sI{2)) 
of the standard loop algebra associated to some fixed divisor of poles |2^ 

D = {{6,,lj),j = 1,...,M,{^,K)}. 



According by ( 5.10| ) the Poisson map 



t,2n 



Cd{sI{2)) is completely defined by the alone 



entry e(A) and by the hamiltonian H related to the linear functional ^x. We put 

_ -"oIKA - Uj) 



e(A) 



^a(^) 



^jk 



K M h 

i=i 



n(A-<5,) 



Res 



A=oo 



(5.16) 



that self-consistent with ( |5.9| ) and (5.13) |14]. Residues and Cjk are the special functions 
of the canonical coordinates {qj}^=i, and variables Uj are curvilinear separated coordinates 
|14|. By using function B residues ej and Cjk are easily restored in canonical variables {(/j}"^^. 
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In the simple poles at A = 5j function e(A) = B"^ ( ^.16 ) has the following residues Cji = 



Parametrization of residues in the higher order poles is discussed in |14]. As an example, we 



consider the polynomial part of e(A) corresponding to a pole at infinity. Let 

K 

B{\) = Y.qK-,\' , (5.17) 



with go = 1 and 



K K K 

e(A)=^e,A^ h{\) =Y,h,\^ , f(X) = Y^f^xj, (5.I8) 

j=0 j=0 j=0 



with ex = ^, = and /i^- = due to the definition of A ( p.9| ) and <I>a (15.161) . Taking 
(|5.13| ) and ( ^.16| ) into account we get 

K-j K-j 

Sj = Qi QK-j-i , = ~ X! QK-j-i , 

i=Q 4=0 

(5.19) 

K-j K-j 

fj — ^ ^ Qi,x QK—j—i,x ^ ^ Qi,xx QK—j—i ; 

i=0 i=0 

where qx = {H, q} and we used the Newton formulae for a product of two sets. Canonically 
conjugate to the coordinates qj momenta pj can be derived from the r-matrix algebra ( p. 10 ) 
|14|. Some first polynomials are equal to 



K = 0, e(A) = 1, /i(A)=0, 

K = l, e(A) = \ + 2qi, -h{X) = pi , 

K = 2, e(A) = + 2Xqi + {2q2 + qf) , 

-h{X) = Xp2 + {pi + P2qi) , , (5.20) 

K = 3, e(A) = X^ + 2X\i + X{2q2 + ql)+2{qs + qiq2), 

-h{X) = X'^P3 + X{p2+P3qi) + iPl+P2qi+P3q2) ■ 



Notice that the kinetic part of the hamiltonian H ( |5.5D has a nondiagonal form in these 
variables 

K 

T = J2PjPK+i-j ■ (5.21) 



Now we turn to the superintegrable systems. For the Taylor projections (M=0) ( |5.7| ) the 
mapping ( |5.6D preserves the property of superintegrability if and only if < i^T according to 
(|5.6| ) and ( |5.16| ). Here AT is a highest order of a pole of the entry e(A) ( ^.161) at infinity and 
is a highest power in projection ( |5.7D . By N > K dynamical r-matrix ( p.8\ ) has the higher 
poles at A = 00 and the corresponding dynamical systems are no longer the superintegrable 
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systems. All these superintegrable systems are related to the special point of r-matrix and 
the associated second Lax matrix A' (5.6) remains a constant in spectral sense under the 
mapping All these superintegrable systems are related to the special stationary flows 

of the KdV hierarchy ( |5.11| ) and genus of the associated spectral curve of L'{X) determined 
by the characteristic equation ( [4.24| ) is no more than number of the degrees of freedom. 

As an example, we present here several superintegrable systems. Let the entry e(A) has 
the simple poles at A = 6j and the K poles at A = infty 



e(A)=Pi,(A)+ 



j=K+l 



X-6, 



where PrW are polynomials given by ( |5.2[1|). For the some first values of K 
corresponding superintegrable potentials in ( ^.51) are 



the 



K = 0, 
K=l, 
K = 2, 
K = 3, 




j=i \ 

j=2 \ «j / 



V = Aqf-8qiq2+J2\^q] + ^j , 

V = -5qf + I2qlq2 - 4gi - Sq^qs + J2 [q] + § 

j=4 \ 



If N = K, then the dressing mapping ( |5.6| ) has the A^ + 1 arbitrary parameters given by the 
function 0(A) = ^jLo^^j^'^ ■ These parameters are related to the canonical shifts of variables 
qj — > qj + Oj and to the common rescaling V — > a^V in the presented potentials. Additional 
integrals of motion may be constructed by the rule (|5.4| ) from the multivariable generating 
function S2(A,/i) (4.28). Of course, for the oscillator all these integrals (5.4) coinside with 
known ones (|3.21| ). 

We may construct the similar superintegrable systems with rational potentials by using 
a general form of the entry e(A) ( ^.161) with the higher order poles and applying more general 
Laurent projection {f).7(). The presented method can be employed to construct superinte- 



grable systems on the other Riemannian spaces of constant curvature |16]. The corresponding 
quantum systems may be obtained by canonical quantization 



Taking into account the r-bracket ( 2.1C| ) one can conclude that the entries e(A) and /(A) 
could play the roles similar to the standard creation and annihilation operators for harmonic 
oscillator 10|. By using the similar transformation of matrices L(X) or L'(X) with matrix 
U ( 3.23 ) at aj = = 1 we can obtain the symmetric representation of these matrices 



e(A) = e(A,ai,a2, . . . ,a„) 



/(A) = e+(A), hiX) = h+iX). 



In this symmetric representation of the Lax matrix the usual method of spectrum-generating 
algebras [§, 19 1 is a part of the standard Bethe ansatz |15, 10|. It should be emphasized, 
that the algebraic Bethe ansatz is a sufficiently universal procedure, which slightly depends 
from particular system in question. It allows us to interpret various concrete models as some 
representations of alone generalized model, which is defined by its r-matrix only. 

As a third example, let us consider the Calogero-Moser systems. It is well known ||2l|, that 
both the Toda models and the Calogero-Moser models are obtained by hamiltonian reduction 
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of the geodesic motion on the cotangent bundle T*G of a Lie group G. For the geodesic 
motion on symmetric spaces of zero curvature the canonical 2-form, the free hamiltonian and 
equations of motion are equal to 



w = tr(d7/ A dx) 
X = y, 



y = 0. 



(5.22) 



For the geodesic motion on symmetric spaces of positive or negative curvature these quantities 
read 



w 



1 . -IX 

-ti{yx ) 



1\2 



(5.23) 



tr(x ^dy A x ^dx) , H 

x = y, y = yx^y. 

The hamiltonians ( p. 22 ) and ( 5.23| ) have the following sets of integrals in the involution 

h = tr(/) and h = tr{yx'^)'' . 

The additional integrals - "projections of angular momentum" (^]^) - are equal to 

Ijk = tr(gp?-i)tr(/) - tr(p')tr(g/-^) . (5.24) 

Here q = x and p = y for the first equations of geodesic motion ( |5.22| ) and g = In x or g = In y 
with p = yx~^ for the second equations of geodesic motion ( [5.23| ). 

In the reduction process the Lax matrices of the reduced system are expressed in terms 
of X by a formula of the type L = zxz~^, where z is some elements in G [21|. For the 
geodesic motion (p. 22 ) associated to the Calogero model with the rational potentials, the 
hamiltonian H ( ^.22]) remains superintegrable and images of integrals ( ^.24 ) are integrals of a 
reduced system ||2^. For the second geodesic motion (5.23) associated to the Calogero-Moser 
model with the trigonometric potentials and to the Toda model, the reduced hamiltonian 
do not commute with the images of logarithmic additional integrals ( |5.24| ). In the quantum 
mechanics whole polynomial algebra of the integrals of motion for the Calogero model has 
been introduced in |17|. 

Let us show as the superintegrable hamiltonian ( ^.22| ) appears in the r-matrix formalism. 
For instance, consider the Euler-Calogero-Moser system P7||. Introduce a set of dynamical 



(5.25) 
(5.26) 



variables {{qj,Pj)}j=i and {fij}ij=i {fij = —fji) together with the Poisson brackets 



{Pj,Qk} = (5. 



{fij^fkl} = 2 i^ilfjk + ^kiflj + Sjkfil + Sljfki) 



In order to have a nondegenerate Poisson bracket it is assumed that the variables fij are 
restricted to a symplectic submanifold of ( |5.26| ). The hamiltonian and the Lax matrix for the 
Euler-Calogero-Moser system |27] are given by 



. N 1 ^ 



{Qi - QjY 



(5.27) 



N ^ { I 

L{\) = Y.Pi^n+ E — 
j=i i,j=i i^j \'^' 



fi j 6 
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with the corresponding r-matrix in the form |^ 



n2(A,/i) 



N 



2 ^ ^jj 




(5.28) 



vy ^ji 



where {eij)ki = hk^ji- In the reduction process this r-matrix inherits the singular point A = 00 



from the initial rational r-matrix. The corresponding superintegrable hamiltonian ( 5.27 ) may 
be defined by (O) with </> = 1/2 • A"^ 



H = ^y 



tiL^{X) 



Res I 



A=oo 



A"^ • trL2(A) 



and the second Lax matrix is independent on spectral parameter 



N 



A = $Atri [rai (A, fi)Li (A)] = ^ _ 

i,j=i ijLj wi Qj) 



The higher flows with (p{X) = 1/k - A ^ ^ in (O) are superintegrable as well [27|. 



6 Conclusions 



We have seen that superintegrable systems closed to geodesic motion can be realizing as 
isospectral flows on coadjoint orbits of loop algebras in framework of the r-matrix formalism. 
All these systems with rational potentials are associated to the special singular point of r- 
matrix. 

Another classical superintegrable systems with an arbitrary number of degrees of freedom 
are the Kepler problem ||l|, |2l| and the rational Calogero-Moser systems . In the proposed 
scheme we can consider a free geodesic motion on the momentum sphere and use the stereo- 
graphic projection with appropriate change of the time variable to study the Kepler problem 
|pO|| . However, this transformation could violate the r-bracket ( ^.10 ) for the corresponding 
Lax matrix. It would be interesting to construct the 2x2 Lax matrix for the Kepler problem 
and for the Kepler-like superintegrable potentials listed in 

This research has been partially supported by RFBR grant 96-0100537. 
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